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Abstract 

The electronic energy band structure of strained and unstrained Si, Ge and SiGe alloys is ex- 
amined in this work using thirty- level k.p analysis. The energy bands are at first obtained with 
ab initio calculations based on the Local-Density-Approximation of Density-Functional Theory, 
including a GW correction and relativistic effects. The so-calculated band structure is then used 
to extract the unknown k.p fitting parameters with a conjugate- gradient optimization procedure. 
In a similar manner, the results of ab initio calculations for strained materials are used to fit the 
unknown deformation potentials that are included in the present k.p Hamiltonian following the 
Pikus and Bir correction scheme. We show that the present k.p model is an efficient numerical 
method, as far as computational time is concerned, that reproduces accurately the overall band 
structure, as well as the bulk effective density of states and the carrier effective masses, for both 
strained and unstrained materials. As an application, the present thirty-level k.p model is used to 
describe the band offsets and the variations of the carrier effective masses in a strained material, 
a Sii-xGe^/Sii-^Gey layer system. 

PACS numbers: 71.15.-m, 71.20.Nr, 71.15.Mb, 71.20.-b 
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I. INTRODUCTION 



The continuous downscaling of MOSFET critical dimensions such as the gate length 
and the gate oxide thickness has been a very successful trend in current manufacturing, as 
testified, e.g., by the ITRS requirements. However, conventional scaling down of MOSFET's 
channel length is declining as the physical and economic limits of such an approach are 
coming closer. Novel solutions are increasingly being used in MOSFET channel engineering. 
The growth of a strained Si layer on a relaxed Sii_ y Ge y buffer layer is a typical technique 
used to improve the electrical performances of MOSFET's device, due to the expected 
enhancements in the carrier mobility 1 ' 2 of such a strained layer. High-performance strained 
Sii-^Ge-j; transistors have been recently obtained and this technology could feature in future 
chip generation with channel size of 32 nm or less 1 ' 3 . When modeling the electrical currents 
of such devices, it is required to take into account the fundamental transport properties of 
the charge carriers that are governed by the structure of the electronic energy bands of the 
strained material. 

Material science computational methods for the calculation of the electronic energy band 
structure fall into two general categories. The first category includes the ab initio methods, 
such as Hartree-Fock or Density Functional Theory (DFT), which calculate the electronic 
structure from first principles, i.e., without the need for empirical parameters. The second 
category consists of far more computationally efficient semiempirical methods such as the 
Empirical Pseudopotential Method (EPM), the Tight Binding (TB) method and the k.p 
method. 

Over the past decades, the Local Density Approximation (LDA) variant of DFT 4 has been 
established as a very powerful tool for studying the elastic properties and the deformation 
potentials of strained semiconductors 5-8 . More recently, the GW many-body-correction to 
the LDA DFT 9 has yielded semiconductors band structures that feature band gap values 
near their experimental values. Ab initio methods are self-consistent methods which utilize 
a variant approach to calculate the ground state energy of a many-body system and thus 
require large computer resources. They can only be used in particular situations of high 
symmetries and are not suitable for calculating the transport properties of large systems 
with confined edge states. 

Unlike ab initio approaches, EPM, TB and the k.p method involve fitting parameters to 
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reproduce the experimental energy band gaps, the dielectric response and the carrier effective 
masses. Over the past three decades, the EPM with spin-orbit (SO) corrections has proven 
to be extremely successful in calculating the electronic band structure of relaxed and strained 
semiconductors with indirect gap 10-14 . Recent works using the TB method have also given 
accurate results 15,16 . In the k.p method, the energy band structure is obtained by a set of 
parameters which represent the energy gaps at T, the momentum matrix elements and the 
strength of the SO coupling. The number of energy bands (or levels) that are effectively 
calculated is related to the precision of the results. Six-level k.p model 17 , eight-level k.p 
model and fourteen- level k.p model 18 describe well the highest valence bands (VBs) and the 
lowest conduction bands (CBs) of semiconductors near the center T of the Brillouin zone, but 
fail to describe the CBs of semiconductors with indirect gap. Low-order k.p Hamiltonians 
need a small number of parameters (typically less than 10), while high-order k.p methods 19 , 
referred by Pollak et al. 20 as "full- zone" k.p methods, require a large number of unknown 
parameters. While it is a straightforward matter to work out the energy-band structure at 
any /c-point in the Brillouin zone once these parameters have been chosen, it requires effort 
and skill to come up with a satisfactory set of parameters 20 . For this reason, the full-zone k.p 
method has been used rarely and for a few bulk semiconductors only, including Si, Ge 19 ' 21 ' 22 
and ct-Sn 20 , while no extensive works were performed for strained semiconductors. 

In this paper, we propose to extend to strained Si, Ge and SiGe alloys the thirty-level 
k.p model that was firstly introduced by Cardona et al. 19 . Cardona's model is based on the 
15 orbital states referred to T. For the first time, the well-known Pikus and Bir 23 correction 
for strained materials has been combined within this thirty-level k.p formalism. The k.p 
parameter optimization strategy is based on a conjugate-gradient procedure that uses ab 
initio simulations but also a large amount of experimental data that is currently available 
for Si and Ge. 

In the first part of this paper, a series of ab initio DFT-LDA simulations that include 
the GW correction and relativistic effects in strained Sii-^Ge^/Sii-^Ge^ systems has been 
performed with a view to complete the experimental data and to establish a reference set of 
energy bands. 

In the second part of this paper, we have determined the k.p model coupling parameters 
and the deformation potentials that fit as close as possible the first principle results, matching 
not only the energy levels and the carrier effective masses, but also the general shape of the 
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band structure of relaxed and strained crystals. A simple interpolation between the k.p 
parameters for Si and those for Ge has been proposed in order to model SiGe alloys. 

In the third part of this paper, a set of comparison is given with experimental data in 
relaxed and strained Si, Ge, and Sii^Ge^ alloys. We show that the present k.p model 
accurately reproduces the overall band structure, as well as the band shifts, the carrier 
effective masses and the Density Of States (DOS) vs. applied strain. A second set of 
comparisons with the widely used Chelikowsly and Cohen non-local EPM 10 has also shown 
a good agreement with k.p simulations. 

II. FIRST PRINCIPLE SIMULATION SETUP 
A. Bulk material 

A series of first principle calculations has been performed in Si and Ge to obtain a 
reference set of energy bands, which can be used later for the k.p model parameters op- 
timization. The first principle data presented in this work have been obtained within the 
LDA variant of the DFT 24 . The present DFT-LDA calculation relies on the pseudopotential 
(PP) approximation, by which the core states are effectively eliminated from the calcu- 
lation. We have used the well-known published Hartwigsen-Goedecker-Hutter relativistic 
separable dual-space Gaussian PP's 25 that uses Cerperley and Alder exchange-correlation 
functional 26 . These PPs include relativistic effects and provide an accurate description of the 
top VB in the near T-region, a critical region with respect to the hole transport properties 
in semiconductors. 

The value of the equilibrium lattice parameter has been calculated by minimizing the 
total energy. Any further LDA calculation has used this theoretical value, instead of the 
experimental one 27 , yielding thus a consistent set of zero-pressure reference data. In Si, we 
found a =5.387 A, and in Ge a =5.585 A. These values agree within 0.75% and 1.33% with 
the experimental values of 5.431 A and 5.658 A 28 , respectively. 

It is known that the band gaps calculated with the LDA method are generally below 
their experimental values. However, the agreement can be greatly improved by the use 
of the Hedin's GW correction 9 . In practice this correction can be applied as a post-DFT 
scheme 29-31 in a non self-consistent way. In the following work, the GqWo correction of bulk 
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Si and Ge were computed on 19 high-symmetry points in the Brillouin zone and added in 
a perturbative manner to the LDA band structure. One remark deserves notice: In spite 
of the G W correction, the theoretical lowest CB typically lie within 0.05-0.2 eV of the 
range of experimental energies observed 29 . In the present work, the GoWo theoretical in- 
direct gaps have been found to be 1.076 eV in Si (located at 84% away from T along the 
T-X direction) and 0.64 eV in Ge, which is underestimated by 8% and 14%, respectively. 
Our results compare favorably with other LDA-GW results found in literature (e.g., see the 
extensive comparison between first principle calculations summarized in Aulbur et al. re- 
view 29 ). Eventhough correctly parametrized, LDA-G W results do not match perfectly with 
experimental data, and significant theoretical work such as vertex correction, self-consistent- 
GW and the exact treatment of exchange term are currently in progress to further improve 
DFT results 29 . These alternative approaches are beyond the scope of the present work. For 
the purpose of obtaining reference set of energy bands that can be used in the development 
of an optimized k.p model, we used the non self-consistent Go Wo approximation to correct 
the band gap problem and we applied a supplementary rigid "squizor" shift of 0.09 eV for 
Si and 0.104 eV for Ge in order to obtain the final reference set of energy bands, the GW 
energy levels 32 . 

Aside from the above studies, the accuracy of the DFT-LDA calculation critically de- 
pends on the manner in which the problem is sampled numerically 33 . We found out good 
LDA convergence (AE k « 0.01 eV) using a basis set of approximatively 1300-1500 plane 
waves, which corresponds to a cut-off energy of around 22 Hartrees. The Brillouin zone was 
sampled on a 6*6*6 four- fold shifted Monkhorst and Pack grid 34 (i.e., 864 Appoints) to obtain 
the charge density. As regard GoWo correction, a satisfactory trade-off between numerical 
convergence and computation time was achieved with a cut-off energy of 8 Hartrees and 
using a large number of bands (> 100) included in the calculation of the self-energy 35 . 

B. SiGe alloys 

In Si^Gex^, where x denotes the relative mole fraction of the two materials, both Si and 
Ge atoms are present in the unit cell. For this reason, we used a 32-atoms tetragonal cell to 
simulate SiGe alloys (a 2-atoms orthorombic unit cell was used for Si and Ge). The Si and 
Ge atoms have been randomly distributed in the super-cell and a structural optimization of 
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the atomic positions in the unit cell has been performed. A linear interpolation between the 
GW correction of Si and Ge was used to correct the band gap. This latter approximation is 
reasonable because the GW correction obtained in Si and Ge are effectively very close 29 . 

The experimental lattice parameter in Si x Gei_ x is well described by the Dismukes's law 36 
according to which a exp = 5.431 + 0.2a; + 0.027a; 2 is a quadratic function of x. A similar 
expression has been obtained from our theoretical results at various x-content: 

a theo = 5.387 + 0.1428a; + 0.0532a; 2 . (1) 



C. Strained material 



Epitaxial Sii^Ge^ layers grown on a relaxed Sii-^Ge^ buffer is studied with a view to 
fitting the deformation potential parameters of the k.p model for strained materials that will 
be presented hereafter. A series of first principle simulations of the electronic band structure 
has been performed for a large range of biaxial strain (up to 4 %) applied perpendicularly 
to the [001], [111] and [110] directions. 

Using continuum elasticity theory, the strain tensor in the layer in the case of [100]-buffer 
writes, 



in the case of [110]-buffer, 
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and in the case of [lll]-buffer: 
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where the longitudinal strain e\\ = a\\/a — 1 depends on the slight difference between the 
cubic lattice parameter in the buffer (a ) and the longitudinal one in the layer (a\\). These 



equations depend on the normal strain e± = —D.e\\. The Poisson ratios D determine the 
displacements of atomic plans along the normal [001], [111] and [110] directions, in order to 
minimize the normal stress. 

The knowledge of the elastic constants of the material is not enough to infer the position 
of the atoms in the unit cell. An additional degree of freedom (that occurs notably for 
shear distortion) must be added to the displacement of the atoms 37 . The internal strain 
parameters £, which measure how the distance between the two atoms in the unit cell 
changes in response to the symmetry-breaking stress, are known experimentally in Si and 
Ge uniaxially strained along the [111] direction 38 (and are quoted in Table I). 

The normal stresses and the inner displacement of atoms in the cell affect the electronic 
and structural properties of the strained crystal 39 . The present calculations, that include a 
structural optimization of the unit cell, are performed as follows: 

(i) a|| is determined from Eq. 1. Eqs. 2, 3 and 4 are used to calculate the shifted Bravais 
lattice R' = (1 + e) • R. 

(ii) The internal strain parameter and the Poisson ratios are calculated by minimizing the 
total energy of the biaxially strained crystal. For searching the equilibrium structure under 
external applied strain en, the total energy is minimized by varying the unit cell length along 
the strain direction and the atomic position in the unit cell. 

(iii) The optimized cell is used later in the DFT-LDA band structure calculations. 

(iv) The GW correction of bulk Si and Ge is used to correct the bandgap problem. This 
choice was motivated by the work of Zhu et al. 40 reporting that there is no quantitative 
difference in Si between the LDA band gap pressure dependencies and the ones from full 
GW calculation. 

The theoretical coefficients D oi, D m and Duo, and the internal strain parameter ex- 
tracted from the previous structural optimizations (step ii) at e\\ — > are reported in Ta- 
ble I. Using continuum elasticity theory, these coefficients write _D oi = 2(Ci 2 /Cn), -Duo = 
(C u + 3Ci 2 - 2C 44 )/(C n + C12 + 2C U ) and D 1U = (2Ci 1 +4Ci 2 -4C 44 )/(Ci 1 +2Ci 2 +4C 44 ), 
where Cn, Ci 2 and C 44 are the elastic constants. Also reported in Table I are the theoretical 
elastic constants obtained from D oi, D m and D 110 together with the above equations. As 
can be seen, the agreement with experimental values is good. 

In Fig. 1, the theoretical Poisson ratios D and the internal strain parameters £ are reported 
as a function of applied biaxial strain e\\. Simulations have been performed in tension in 
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Si (right) and in compression in Ge (left) for the three previously mentioned growth cases. 
As can be seen, for the [001]-growth case D does not depend on e\\, while for the [111] and 
[110]-growth cases both D and £ change (up to ~ 15%) with applied strain. The former 
cases is consistent with results of Ref. 6 concerning strained Ge layers grown on [001]-cubic 
Si. Changes in £ with applied strain have also been reported in Refs. 5 ' 8 ' 41 as well as elastic 
constant strain-dependency 39 . As shown later, the energy band shifts notably depend on the 
deformation applied on the crystal and thus the results shown in Fig. 1 have been carefully 
accounted for. For instance, we found out that the VB energy shifts in Si grown on a [111]- 
Ge buffer are overestimated (by roughly 10%) when both £ and D are kept constant (this 
behavior is inversed in Ge grown on [11 1]-Si buffer). 

III. DEVELOPMENT OF AN OPTIMIZED k.p MODEL 
A. k.p parameters for bulk materials 

The k.p formalism using Zinc-Blende T-centered Bloch function basis U[ k (r) = 
^ n C l n u n o(r) leads to the secular k.p equation of the undeformed crystal 17 : 

n 

+ S-<«n'o|p|«nO»C^ = 0, (5) 

where E® are the eigenvalues at T. The fifteen T-states from group O h determined by 
Cardona et al. 19 (shown in Table II) is our starting point. The number of independent non- 
diagonal matrix elements in Eq. (5) can be reduced to ten in Si and Ge using group-theory 
selection rules 19 ' 46 . The SO coupling terms are introduced in the usual way 17 ' 20 leading to 
a 30x30 k.p matrix (See Appendix A). One notes incidentally, that since this approach 
does not use renormalized Luttinger-like parameters, it is different from other lower-order 
k.p models such as the fourteen-level k.p model of Ref. 18 and the twenty-level model of 
Ref. 47. 

As written in appendix A, the thirty-level k.p model depends on seven T-centered eigen- 
values, four SO coupling coefficients and ten matrix elements. In the development of the 
present optimized k.p model, we attempted to fit experimental electronic properties of Si, 
Ge and SiGe alloys as close as possible. Currently, there is insufficient detailed experimen- 
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tal information about the band energies to accurately determine all the k.p parameters, 
particularly at high energy (> 5 eV). For this reason, we adopted a mixed approach us- 
ing experimental data when available and ab initio results otherwise. These coefficients 
were fitted using a conjugate-gradient procedure. Satisfactory convergence was determined 
through a least-square error function between k.p eigenvalues and GW results evaluated 
on a dense 1000) set of Appoints in the Brillouin zone. A particular care has been paid 
in the near-r region and at the CB minima in order to obtain accurate description of the 
curvature masses and Luttinger parameters. We also tried as far as possible to reduce the 
discontinuity in energy at the K, U equivalent points 20 . Due to missing high energy (220) 
bands in the present k.p eigenvalues at K and U differ from several meV (~7 meV). The 
coupling between (220) states to the other thirty-lowest energy states are naturally more 
pronounced near K, where the lowest (220) eigenvalues are as low as ^ 4.5 eV. For com- 
parison, the lowest (220) eigenvalues are > 26 eV at T, > 12.5 eV at X and > 12 eV at L. 
After several attempts to remove this discontinuity by changing the k.p parameters values 
only, we found out that no satisfactory trade-off between accuracy and continuity in K, U 
could be obtained. For this reason, the present k.p model has a discontinuity in K, U 48 (as 
in Ref. 19). 

Following Pollak et al. 20 , the SO strength between and r 25 ; states were determined by 
imposing the highest VB to be degenerated at X. A 25 i = 44 meV for Si and A 25 ; = 290 meV 
for Ge are known by experiments, while A 25 « and A 15 were obtained from first principle 
simulations (see Table II). 

The T-centered eigenvalues and the coupling parameters obtained from our procedure are 
listed in Table II and Table III. One should mentioned that there are two sets of published 
k.p parameters in Si 21 ' 22 and one set in Ge 21 based on the early work of Cardona and 
Pollak 19 . Although Cardona and Pollak parameters set 19 provides an accurate description 
of the main CB minima and the top of the VB, it has the limitation of not including the 
SO coupling (due to computational limitations in the mid-sixties). This was recently done 
in Si and Ge 21 together with a new set of parameters. Unfortunately in Ge, the proposed 
set of parameters failed to improve Cardona and Pollak 19 one in so far as none of the VB 
reach the edge of the Brillouin zone with zero slope (or average slope) as required by crystal 
symmetry and the L- valleys minima are very distant from the Brillouin zone edge. This later 
severe drawback makes this model unappropriate for application to transport properties in 
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nanostructure. This was not the case with Cardona and Pollak 19 parameter set, but also with 
the present one. In comparison to this former set of parameters, our optimization strategy 
based on ab initio reference set of energy bands brings additional informations (notably at 
high energies) and slightly improves the accuracy of carrier group velocity at the first and 
second CB minima, but also certain energy gap values in both Si and Ge. Parameters listed 
in Table II and Table III are slightly different from Cardona and Pollak 19 ones. The main 
differences 49 come from the fact that the T— eigenvalues used in the present k.p model 
differ at high energy (E > 3.5 eV), and that non-local effects have been accounted for in the 
present model 50 . 

The Si and Ge BS obtained using the present k.p model are compared to GW first 
principle simulations in Fig. 2 and Fig. 3. The overall quality of the fit is excellent: The 
difference in band energies between our semiempirical values and those used for the fit was 
typically less than 0.01 eV for the principal band gaps, and under 0.3 eV at other high 
symmetry points. 

Further comparisons are shown in Fig. 2 and Fig. 3 with the widely used Chelikowsky 
and Cohen non-local EPM 10 . As can be seen, the straightforward application of the EPM 
yields a BS that is in very good agreement (up to more than 6eV) with the results of the 
much more complex GW ab initio calculation. This excellent agreement means that both 
methods are consistent and leads further support to the GW reference set of energy bands 
used for the k.p parameters optimization. 

B. SiGe alloys parameters 

The virtual crystal approximation was used to extend the k.p results to SiGe alloys. 
A quadratic interpolation between Si and Ge parameters is proposed in Tables II and III. 
Due to centro-symmetry breaking SiGe alloys do not belong to group O h . A supplementary 
SO coupling term 51 and two purely imaginary coupling terms from group T d52 (not shown 
in appendix A) have been introduced in Table III. These interpolation coefficients were 
determined in order to account for the GW Sii^Ge^ band structures for various re-content 
values. Fig. 4 shows the Sii_ x Ge x SO splittings and band gaps as a function of x calculated 
using the present k.p model, but also with EPM and GW model. The present k.p model 
predicts a crossing between A-valley and L-valley minima at x=0.84 which is consistent 
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with the experimental data 1 . With these new parameters and interpolation coefficients, one 
obtains a good agreement between the k.p model and GW results as it is testified by the 
Si .5Ge .5 band structure shown in Fig. 5. 

C. Strained materials: deformation potentials 

The analysis of the structure of strained semiconductors using k.p Hamiltonian has been 
initially proposed by Pikus and Bir 23 . Authors derived the first-order k.p perturbation 
terms arising from straining the semiconductor in question. This general expression has 
been widely used in six- level k.p Hamiltonian to analyze the effect of deformation on the 
hole energy spectrum. For the first time, we applied Pikus and Bir formalism in the case of 
the thirty- level k.p Hamiltonian. The perturbation term to be added to Eq. 5 writes 23 : 

W n > n = -— J2i j e ijki- < U n 'o\Pj\u n0 > 
171 ,J 

+€ij < w n / |Hij|w„o >, (6) 

where i,j stand for x,y,z. The first term of Eq. 6, accounts for the interaction between 
the strain and the momentum of the carriers. The deformation potential operator = 
— + ^i( r ) describes the change in the potential and the kinetic energy of carriers due to 
the strain itself. The unknown non- vanishing deformation potentials at Y < ti n 'o|— ij|w n o >, 
listed in Table IV, were determined from group theory selection rules taking into account 
the piPj contribution only. This choice is motivated by the fact that in the deformed-ion 
approximation the Vjj term with rhombic or tetragonal symmetries vanishes 23 . 

The perturbation matrix W to be added to the thirty k.p matrix is shown in appendix 
B. Accurate knowledge of the deformation potentials at T for all the thirty lowest energy 
bands is required for the construction of the k.p model in strained material. However, 
only deformation potential for the T 2 5i states are known experimentally 6 . The deformation 
potentials, needed in Eq. 6 and listed in Table V have been fitted using a procedure similar 
to the one used in bulk materials. We have used a least square optimization procedure on 
GW energy levels calculated for various strain tensors, including shear distortions. Special 
attention has been given to the respect of time-reversal symmetry degeneracy at Brillouin 
zone edge. In addition to space-group symmetry operations, the Hamiltonian of an isolated 
centrosymetric crystal exhibits time-reversal symmetry. Therefore, additional degeneracy 
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among eigenvalues may be determined using the Kramer theorem and the Wigner rule. 
Using that rule, Ma et al. 15 have obtained the additional informations about degeneracy in 
Si and Ge for the [001], [111], and [110] growth cases that have been accounted for in the 
present k.p model. 

Typical results of our fit procedure are presented in Fig. 6 (and Fig. 7) for biaxially 
strained Si (Ge) layers on [001]-oriented cubic Ge buffer (and Si buffer, respectively). Further 
results are shown in Fig. 8 (and Fig. 9) for strained Si (Ge) layers on [lll]-oriented cubic Ge 
buffer (and Si buffer). The band structures calculated along various directions in reciprocal 
space using the present k.p model are compared to EPM 53 and GW results. 

Under uniaxial strain along [100], the minima of the L- valleys remain equivalent. How- 
ever, due to crystal symmetry lowering by biaxial-strain, the VB degenerate levels are split 
and equi-energy lowest CB A-valleys are split into four A 4 and two A 2 - valleys. On the con- 
trary, uniaxial strain along [111] leaves the A-valleys equivalent whereas valleys at second 
minima are split (two Z-valleys along [111] and six L-valleys along [111], [111] and [111]). 

Another important effect of the strain is observable in the near-r region: In the case of 
the strained Si (tensile-biaxial-strain), the heavy- hole band "crosses" the SO band, while 
in the case of strained Ge it is the contrary (compressive-biaxial strain); The SO-hole band 
"crosses" the heavy-hole one. 

IV. RESULTS AND COMPARISON TO EXPERIMENTAL DATA 
A. Energy band gaps 

We address now the question, how well do the energy levels (vide post for the effective 
masses) calculated using k.p agree with the experiments and theoretical results. In Tables II 
and VI comparison is given for the cases where experimental data seem to be well established. 
To this end, we have used the large set of experimental values summarized by Chelikowsky 
and Cohen 10 , Landolt-Bornstein 28 , Hybertsen and Louie 45 , Aulbur et al. 29 , but also the 
most recent optical measurements in Si 43 and Ge 44 . One notes that this assignment is in 
some cases somewhat tentative. The interpretation of the experimental peaks and critical 
points is hindered when energetically close transitions take place, e.g., due to SO splitted 
bands. Moreover, the critical points might originate from transitions close to, but not 
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necessarily exactly at high symmetry points in the Brillouin zone. Aside from recent inverse- 
photoemission and photoemission experiments which have addressed certain high energy 
bands in Si and Ge 42 , the complex shape of the different energy bands cannot be easily 
measured directly on its whole structure. 

It is clear from Tables II and VI that the overall agreement between theoretical and 
experimental band gaps is good. In particular, the present k.p model predicts indirect gaps 
of 1.17 eV in Si and 0.747 eV in Ge. 

B. Curvature masses, Luttinger parameters and DOS 

The A-electron and L-electron effective masses were obtained from the second derivative 
of the CB energy with respect to wave vector along various directions away from the valleys 
minima. For the VB, the bands are extremely non-parabolic, and so the effective masses 
could not be evaluated by the method described above. Instead, we used the six-level k.p 
Dresselhauss-Kip-Kittel model 17 that depends on three Luttinger parameters, the values 
of which have been fitted using a conjugate- gradient optimization. This optimization is 
based on a least square error between the curvature masses along the [001], [111] and [110] 
directions obtained with this 6-level k.p model and with EPM, k.p or GW models. The 
present extractions technique of the curvature masses and Luttinger parameters, although 
different to the one proposed in Ref. 14, gives similar results (e.g., the Chelikowsky and 
Cohen-based EPM values obtained in Ge 7x=9.563 , 7 2 =2.77, 7 3 =3.91 can be compared to 
our results listed in Table VII). The theoretical effective masses and Luttinger parameters 
are listed in Table VII and compared to experimental data. All methods provide reasonably 
good agreement with experimental values. One should mentioned nevertheless that the 
present k.p model as well as Chelikowsky-based non-local EPM give rather disappointing 
results for the VB Luttinger parameters in Ge. Even if the accuracy of the Luttinger 
parameters is improved with the present k.p model, it still underestimates their values by 
about 20%. When one expresses the Luttinger parameters in terms of matrix elements 17 , 
it becomes clear why it so: They contain the term Pj Ep ■> the value of which is relatively 
small in the present work 50 . 

Accurate description of the DOS are key features for accurate carrier density modeling and 
realistic transport models in semiconductors 58 ' 59 . They are also good checks for the quality 
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of the present k.p model. Indeed, the DOS not only depends on the band energies but also 
on their gradient with respect to wave vector (group velocity). The DOS is obtained using 
the Gilat and Raubenheimer procedure 60 . We applied exactly the same algorithm using the 
k.p (solid lines) the EPM (dashed lines) and the GW models (dotted lines), respectively. 
As can be seen in Fig. 10, the agreement between models is excellent for the VB but also 
for the CB from -5 eV up to 5 eV. 

C. Energy band shifts in strained Sii -^Ge^/Sii - y Ge y systems 

The calculated energy shifts of the main VB and CB extrema in Si, Ge and Sii„ y Ge<, are 
shown in Fig. 11, 12 and 13 for material grown on [001], [111] and [110] Sii^Ge^, buffers 
(the energy scale has been fixed by setting arbitrarily to zero the top of the VBs). In 
these figures, comparison is shown between first principal simulations, EPM, the present 
k.p models and the close agreement between methods can be noticed. In particular, the 
A- valley and L- valley band shifts are correctly modeled with the present k.p analysis. 

An important effect of strain can be seen for shear distortions ([111] and [110] growth 
cases) in Figs. 11. Several T— eigenvalues shifts exhibit a non-linear relation vs. y. This can 
typically be inferred from symmetry lowering due to strain. Indeed, several /c-independent 
coupling terms (Wr terms in the W30X30 perturbation matrix) occur in strained materials 
(e.g., between the and the T 2 n states originally non-coupled in the relaxed k.p Hamil- 
tonian). 

From the energy-band shifts, we have calculated the deformation potentials which reflect 
the variation of an individual band energy as a function of applied strain. Following the 
notation of C. G. Van Walle 7 , the experimental values found in strained Si and Ge for the 
splitting of the top of the VBs and for the lowest CB minima shifts are quoted in Table V. 
The deformation potentials a v , b v and d v refer to the hydrostatic, the splitting and the 
shear deformation potentials of the VBs at T, respectively. Similar expressions apply for the 
CB, but only a c is experimentally addressed so far in Si and Ge 61 . Theoretical deformation 
potentials are consistent with experimental ones excepted for d v in Ge that is overestimated 
by ~ 8%. This result can be inferred from the slight overestimation of the theoretical internal 
parameter £ shown in Table I. 

Besides deformation potentials at T, other deformation potentials have been experimen- 
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tally determined at the lowest CB minima, typically along the L-direction for Ge and X- 
direction for Si. In Table V, S^, S£ are the hydrostatic CB deformation potentials, while 
S£, are the splitting deformation potential. In general, individual quantities are difficult 
to measure because they are referred to an absolute scale. The case of the VB hydrostatic 
deformation potentials a„ is worth mentioning: Accurate value is hard to obtain from ex- 
periments or from theory. In this work, we did not calculate a„, we fixed its value following 
Fischetti and Laux results 13 : a v =2 for Si, a v —2.1 for Ge and a v = 2 + O.lx for Sii_ x Ge x . 

D. Carrier masses in strained Sii_. E Ge. E /Sii_.yGe, / systems 

Changes in effective masses in case of biaxially strained Si and Ge have been been reported 
by Rieger and Vogl 14 , Gell 11 and by Fischetti and Laux 13 using EPM. In this section, we 
generalize EPM results to shear distorted crystals and to k.p and GW methods. This is to 
our knowledge the first ab initio calculations of the effective masses vs. strain reported so 
far in strained Sii_ :r Ge ;E /Sii„.yGe y systems. 

Figures 14 and 15 show the strained Si and Ge electron curvature masses as a function of 
y-content in the Sii-yGe^ buffer for various orientations. According to previously reported 
EPM results 11 ' 13 ' 14 , electron curvature masses change with applied strain (up to 100%). 
Simulations performed using semiempirical methods (k.p EPM) are consistent with first 
principal results, even though k.p tends generally to underestimate the changes in masses 
vs. strain. 

For biaxially strained Si, the longitudinal and transverse curvature masses depend on 
the A-valley (Fig. 14). For the [111] and the [110] growth cases, the curvature masses 
significantly increase with applied strain. The effects of strain on the L-valley curvature 
masses are generally less pronounced than for the A-valley, although variations up to 20% can 
be observed. Similar conclusions apply for the A- valleys and the L- valleys of Ge (Fig. 15). 

We found out that the A-valley minima positions in reciprocal space also change with 
applied strain. This behavior is generally more pronounced when shear distortions are 
applied. For instance, the A-valley minimum distance along T-X direction changes from 
84% in bulk Si up to 97% in strained Si on [lll]-oriented Ge buffer. This latter point and 
the changes in shape can be seen in the three-dimensional (3D) surface plot at thermal 
energy (| kT) shown in Fig. 16. The A-valleys in the first and in the second Brillouin 
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zones are shown respectively along the T-X, T-Y and T-Z directions for bulk Si and for the 
[100] and the [lll]-growth cases. The A-valleys 3D-surface in the first and in the second 
Brillouin zones are clearly separated in the bulk and in the biaxially strained crystal, while 
they "merge" into a single 3D-surface for the [lll]-strained crystal. 

For the reason mentioned previously, the VB curvature masses cannot be easily evalu- 
ated from parabolic fit. Instead, we calculated the DOS effective masses at thermal energy. 
Figures 17 and 18 show the Si and Ge hole DOS effective masses at 300 K as a function 
of y-content in the buffer. In the unstrained crystal, the heavy-hole DOS mass is approxi- 
mately three times larger than the light-hole mass. This difference decreases as soon as the 
degeneracy at T between the heavy and the light hole is removed. One notes that the results 
shown in Figs. 14, 15, 17 and 18 are consistent with Fischetti's EPM results 13 . 

V. DISCUSSION 

We have presented a thirty-level k.p model for bulk and strained SiGe, the parameters of 
which have been optimized using a conjugate- gradient procedure on a reference set of energy 
bands obtained with first principle calculations. The first principle simulations have been 
validated through an accurate comparison with experimental published results. For bulk 
SiGe, a set of comparisons with experimental data has shown good agreements for the main 
band gaps values (Table II and Table VI), the carrier effective masses and Luttinger param- 
eters (Table VII). For strained materials, we have benchmarked our first principle results 
using the deformation potential theory applied to the specific case of Si and Ge epitax- 
ial layers grown on [001], [111] and [110]-oriented relaxed buffers. Theoretical deformation 
potentials (Table V) have been found to be consistent with experimental ones, within the 
experimental error. A second set of comparisons has been performed using the widely used 
Chelikowsky and Cohen non-local EPM 10 , including relativistic corrections. 

As it is presented in this paper, we obtained good agreements between the present k.p 
model and first principle simulations. Since the k.p parameters have been fitted on first 
principle data, one might describe these excellent results as a puppy chasing its tail. This 
is to some degree true, but our procedure yields additional pieces of information: 

(i) Using first principle simulations, important quantities not addressed experimentally, 
such as high energy levels and effective masses at the second CB minima have been taken 
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into account in the present k.p model. 

(ii) SiGe compounds have been modeled using interpolation functions between Si and Ge 
k.p parameters. Additional SO splittings and energy degeneracy removal (due to centrosym- 
metry breaking) have been accounted for with supplementary coupling k.p parameters fitted 
on first principle data. 

(iii) State of the art "full-band" Monte-Carlo simulators for transport properties in 
strained Si/Ge devices are notably based on EPM 13 ' 58 ' 59 . Because simulations with the 
present k.p model basically give the same results than non-local EPM, the k.p model can 
be used as is into a EPM-based MC simulator. It should be noted that the computational 
burden of the k.p calculations is impressively reduced in comparison to EPM simulations: 
The CPU time to obtain a complete band structure is approximatively two orders of mag- 
nitude less than the CPU time needed by EPM 62 . Further comparisons along that line are 
presented elsewhere 59 . 

(iv) Concerning the effective electron masses, our calculations have confirmed Fischetti 
and Laux non-local EPM results 13,14 which have shown that the curvature masses appear 
to be strongly dependent on the strain. Although never measured directly, recent transport 
simulations in inversion layers using "full-band" MC 63 have suggested that effective mass 
change due to strain should be involved in the mobility variation in strained Si and Ge. For 
the effective masses at the second CB minima in energy (i.e., L- valley in Si and A-valley in 
Ge), similar good agreements between EPM and GW results have been found. It should be 
noted that these valleys could play a significant role in the transport properties of strained 
semiconductors, particularly for Ge, in which A-valley and L-valley are separated by only 
200 meV. It has been shown in Section III, that biaxially strained Ge on Sii-^Ge^ buffer 
exhibits a cross-over between A-valley and L-valley minima for y < 0.5. 

For the reasons mentioned above, the present highly optimized k.p parameters set im- 
proves the published k.p parameters set of Refs. 19 and 21 for bulk materials and extends its 
predictions to SiGe alloy. Furthermore, the Pikus and Bir perturbative treatment of strain 
was for the first time evaluated in the Cardona and Pollak thirty- level k.p formalism 19 . We 
have shown notably that this correction captures the main feature of strained-crystal band 
structures, such as energy shifts versus strain and effective masses change due to strain. 
One should mention that this is not the case with the recently published twenty-level k.p 
model for biaxially strained Si and Ge 64 , in which important contributions to the strain 
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perturbation matrix have been omitted. We found out that the behavior of the CBs in 
strained semiconductors strongly depend on the and Wr,r terms in the perturbation 
matrix. Neglecting these terms, as in Ref. 64, leads to a large underestimation of the CB 
equienergy valley splitting and effective masses changes versus strain. Moreover, VB time 
reversal symmetry at X and correct band-shifts at L- valleys versus strain cannot be obtained 
without these contributions. 

VI. CONCLUSION 

In this paper, we have developed a highly optimized thirty-level k.p model for strained 
Si, Ge and SiGe alloys. 

A series of ab initio DFT-LDA simulations that include GW correction and relativistic 
effects in Si, Ge ans SiGe alloys has been performed with a view to obtaining informations 
not addressed by experiments. Once a reference set of energy bands has been obtained, we 
have optimized the k.p model parameters using a conjugate- gradient procedure in order to 
fit as close as possible first principle results, but also carrier effective masses and luttinger 
parameters. A simple interpolation between Si and Ge k.p parameters has been proposed 
in order to model SiGe alloy. 

The electronic structure of strained Sii-zGe^/Sii-yGe,, systems has been studied using 
first principle simulations. For the first time, the well-known Pikus and Bir correction 23 for 
strained materials has been examined within this thirty- level k.p formalism. The deforma- 
tion potentials have been obtained from first principle simulations in order to fit the shifts 
of the thirty lowest energy levels at T vs. applied strain, but also the general shape of the 
band structure of the strained crystal. 

Finally, the present k.p model has been validated through an accurate set of comparisons 
with experimental data in relaxed and strained Si, Ge, and Sii^Ge^ alloys. A second set 
of comparisons with first-principal simulations, but also with the widely used Chelikowsly 
and Cohen non-local EPM 10 has also shown a good agreement. The present k.p description 
of strained Si, Ge and SiGe accurately reproduces the overall band structure, as well as the 
band shifts, the carrier effective masses and the DOS vs. applied strain. 
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APPENDIX A: k.p MATRIX FOR RELAXED MATERIALS 

The thirty-level k.p matrix (Eq. 5) for relaxed materials writes: 



tt30 
H k.p 



H 



2x2 



pill jp2x& 

















pn H 2x6 


jpQxQ 


R'H 6 k x4 








Q'tff 6 


pi H 6x2 


m° r 

1 25'™ i 1 25' ; 




jj4x4 
r i2' 




















jj2x2 

r i u 





Tiff 6 














tt2x2 


T'if 2x6 

rr6x6 

Tl5 





jp2x2 
F 2" 




p H 2xS 
jjQxQ 

r 25" 



(Al) 



whose diagonal blocks read: 



jj6x6 
jj4x4 



2x2 



dtag (Er + ^) + # r S 



(A2) 



where k 2 = k 2 + k 2 + k 2 , and diagi) stands for the diagonal matrix. E r is the eigenvalue of 
the state labeled by T, as listed in Table II. The coupling constants (P,P',P"', etc ...) are 
listed in Table III. is the SO matrix, which depends on the SO coupling parameters 

listed in Table III: 

-1 -i 1 
% -1 -% 
—1—1 i 
—1—1 i 

-% -% -1 

1 i -1 

The non-zero k.p blocks write: 



H 



so 



Ar 
3 



(A3) 
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Jj6x6 



jj4x6 



Jj2x6 



k z k y 

k z k x 
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k z k x 

k v k r 



— A/3fc 2 

— k z 









\piky —\f?>k z 

2A/^> ^y ^ z 





k x ky k z 



(A4) 



We have finally included a SO coupling term between the r 2 5^ and the T 2 ^i states and 
fitted the coupling strength A r r in order to respect the time reversal degeneracy at X. 



APPENDIX B: k.p MATRIX FOR STRAINED MATERIALS 

The perturbation matrix for strained materials (Eq. 6) writes: 



2x2 

- 2'" 
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1 25'" 
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"25'' 



There are two types of coupling terms in the matrix described by Eq. Bl; k-independent 
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terms (labeled Wr) coming from the second term in Eq. 6 and terms that are linear in k 
(labeled Wk) coming from the first term in Eq. 6. 
The k-independent W T blocks write: 



6x6 



^3X3 
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w^ 6 = h 
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(B2) 

The deformation potentials (/, n, m...etc) are listed in Table IV (the coefficients not men- 
tioned in the Table are set to zero). Group theory considerations allow to write the five 
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coefficients A, B,C, D, E as a linear combination of four coefficients: 



A = 


6 (612 - d 12 ) 


B = 


3 (ai2 + 612 - 2ci 2 ) 


C = 


2 (2oi 2 - 4ci2 + 612 + di 2 ) 


D = 


56i2 - 2ci2 - Adi2 + 012 


E = 


\/3 (2ci2 - 2c?i2 - a i2 + 612 



The first term in Eq. 6 gives rise to an additional non-diagonal k-independent coupling 
between states of the same polarity. These blocks and the corresponding deformation poten- 
tials have been labeled using a double subscript notation (e.g., Wr 12 ,r 2 „)- For simplicity, we 
dropped the double subscript notation for the coupling between identical states (e.g., Wr 12 )- 

Finally, the k.p coupling terms due to the second term in Eq. 6 write: 
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(B4) 



where i stands for x, y and z. 
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FIG. 1: Theoretical internal strain parameter £ and theoretical normal strain coefficient D as a 
function of longitudinal biaxial strain applied perpendicularly to the normal [001], [111] and [110] 
directions. The coefficients D are normalized to the experimental values listed in Table I. 
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FIG. 2: Bulk Si electronic band structure obtained using thirty-level k.p model, EPM and GW 
calculation. 




FIG. 3: Bulk Ge electronic band structure obtained using thirty-level k.p model, EPM and GW 
calculation. 
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FIG. 4: Relaxed Sii_ x Ge x band gaps and SO splittings obtained using thirty-level k.p model, 
EPM and GW calculation. 
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FIG. 5: Bulk Sio.sGeo.s electronic band structure obtained using thirty-level k.p model and GW 
calculation. 
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FIG. 6: Electronic band structure of strained Si layer grown on [001]-oriented cubic Ge buffer. 
Simulations with thirty- level k.p model (solid lines), EPM (dashed lines) and GW (dotted lines). 
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FIG. 7: Electronic band structure of strained Ge layer grown on [001]-oriented cubic Si buffer. 
Simulations with thirty- level k.p model (solid lines), EPM (dashed lines) and GW (dotted lines). 
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FIG. 8: Electronic band structure of strained Si layer grown on [lll]-oriented cubic Ge buffer. 
Simulations with thirty- level k.p model (solid lines), EPM (dashed lines) and GW (dotted lines). 
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FIG. 9: Electronic band structure of strained Ge layer grown on [lll]-oriented cubic Si buffer. 
Simulations with thirty- level k.p model (solid lines), EPM (dashed lines) and GW (dotted lines). 
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ENERGY (eV) 

FIG. 10: The DOS of Si, Ge and Sio.sGeo.s alloy: Comparison between thirty-level k.p (solid lines), 
EPM (dashed lines) and GW (dotted lines). 
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FIG. 11: Calculated VB and CB shifts of a strained Si layer as a function of y-content in the 
Sii-yGe^ buffer: k.p (solid lines), EPM (dashed lines) and GW (dotted lines) simulations per- 
formed for various buffer orientations. 
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FIG. 12: Calculated VB and CB shifts of a strained Ge layer as a function of y-content in the 
Sii-^Ge^ buffer: k.p (solid lines), EPM (dashed lines) and GW (dotted lines) simulations per- 
formed for various buffer orientations. 
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FIG. 13: Calculated VB and CB shifts of a strained Sii^Ge^ layer grown on a Si buffer: k.p (solid 
lines) and GW (symbols) simulations performed for various buffer orientations. 
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FIG. 14: Strained Si electron longitudinal and transverse curvature masses as a function of y- 
content in the Sii-yGe^ buffer: k.p (solid lines), EPM (dashed lines) and GW (dotted lines) 
simulations performed for various buffer orientations. 
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FIG. 15: Strained Ge electron longitudinal and transverse curvature masses as a function of y- 
content in the Sii-yGe^ buffer: k.p (solid lines), EPM (dashed lines) and GW (dotted lines) 
simulations performed for various buffer orientations. 
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a) Bulk Si 



b) Si on [001]-Ge 



c) Si on [11 1]-Ge 



FIG. 16: 3D-surface plot at thermal energy of the lowest CB across the Brillouin zone edges 
along the [100], [010] and [001] directions. Simulations performed for a) bulk Si, b) strained Si on 
[001]-oriented cubic Ge buffer, and c) strained Si on [111] -oriented Ge buffer. 
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FIG. 17: Strained Si hole DOS effective masses as a function of y-content in Sii-^Ge^ buffer: k.p 
(solid lines) and GW (symbols); T=300 K. 
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FIG. 18: Strained Ge hole DOS effective masses as a function of y-content in Sii-yGe,, buffer: k.p 
(solid lines) and GW (symbols); T=300 K. 
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TABLE I: Elastic coefficients and internal displacement. 







Si 


Ge 






Exp. 


LDA" 


Exp. 


LDA a 


Dooi (GPa) 


0.776 


0.795 


0.7513 


0.711 


Olio (GPa) 


0.515 


0.527 


0.4498 


0.42 


Dm (GPa) 


0.444 


0.461 


0.3711 


0.343 


Cn (GPa) 


167.5 d 


168. 3 6 


131.5 d 


132.8 6 


C12 (GPa) 


65 d 


66. 8 b 


49.4 d 


46. 8 6 


C44 (GPa) 


80.1 d 


79.9 b 


68.4 d 


66.57 6 




0.54 c 


0.536 


0.54 c 


0.495 


"present work; 


b clastic constant calculated from the D values; 


C C. S. G. Cousins ct al., J. Phys. 


C 20, 29 



(1987); d H. J. McSkimin, J. of Appl. Phys. 24, 988 (1953), cited by Rcf. 7. 
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TABLE II: Eigenvalues and SO splittings of T-centered states. Symbols of Ref. 19. The r 25 ,; 
eigenvalue is arbitrarily set to zero. A symbols refer to the state's SO splittings. All energies are 



expressed in eV. 



States 


Si 






Ge 






at r 


Exp. 


EPM a 


G W a Exp. 


EPM a 


G W a 


k.p e 


rn 


-12.4+0.6 6 ; -11. 2 i ; 


-12.36 


-11.489 -12.6+0.3 6 ; 


-12.624 


-12.638 


-12.7-0.18a; 




-12.5±0.6 C ; -11.4' 




-12.9±0.2 h 








A 25" 


0.044 c 


0.044 


0.0499 0.296 c 


0.297 


0.312 0.044+0.22+0. 052:c 2 


Fib 


3.4 C ; 3.35-f; 


3.406 


3.204 3.006*; 3.206*; 


3.279 


3.1 


3.335-0.2222 




3.05 d 




3.16 9 ; 3.25 c 








Ais 


0.04 c 


0.037 


0.037 0.200 c 


0.205 


0.227 


0.033 + 0.157a: 


r 2 , ; 


4.15'; 4.1 d ; 


4.062 


3.96 0.89' 1 ; 0.90* 


0.861 


0.715 


4.15-3.262 




4.185 c ; 4.21 1 












r i u 




7.561 


8.308 


6.072 


6.82 


8.4-1. 62 


ria' 




9.371 


8.451 


8.665 


9.925 


8.54+1.762 


r 2B /« 




12.203 


11.41 


11.334 


11.193 


11.7-0.34x 


A 25 ,u 




0.009 


0.012 


0.0558 


0.029 


0.012+0.032 


Y 2 ,u 




13.3 


15.41 


12.97 


14.086 


15.8-1.82 


"present calculations; 6 as 


presented 


in Ref. 10; c as presented in Ref. 28; d Rcf. 42 


; e as used in the present 



k.p model; -^Ref. 43; 9 Ref. 44; ''as presented in Ref. 45; *as presented in Ref. 29. 
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m A T~i T T — 1 III at i , fill* i / \ 1 

lABLb 111: Matrix elements or the linear momentum p (a.u.J used 


in the present k.p model. 


Symbols of group O h are taken from Ref. 19. Other symbols for Sii-^ 


Ge x (0 < x < 1) belong to 


group T d . 




Matrix elements (a.u.) 


Sii— x Ge x 


p = ^2 5 n\p\r 2 n) 


1.22-0.034jj 


Q= ^< r 25"lp|ri 5 > 


1.0679+0.0068a; 


R= ^<r 25 ^|p|ri2'> 


0.5427+0.0884:r 




0.156-0.00812: 


*"=£<r 2B ,»|p|r 2 „) 


-0.008+0.078x-0.05x 2 


Q'= ^<r 25 ,»| P |r 1B ) 


-0.6555-0.1052o: 


R'= £<r 25 ,«| P |r ia ,) 


0.8342-0.0126x 


P"'=^<r 25 ,u|p|r 2 ,„> 


1.425-0.02633: 


r= A (ri „| p |r 1B ) 


1.166-0.0247i>0.04:r 2 


t'= A( ri! |p|ri 5 > 


0.29+0. 08a: 


s= ^(r 15 | P |r 2/; > 


-i0.1a:(l - x) 


s'= A (ri5 |p|r 2 ,.> 


i0.3a;(l - x) 


SO coupling strength (eV) 


Sii-^Gea; 




0.022+0. 198x 


A r 15 ,r 25 , i 


0.04jj-0.04x 2 
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TABLE IV: Strain perturbation matrix coefficients expressed in eV. 



symbols 




symbols 




symbols 


oil — jUj 


symbols 


Oil — x^^x 


25 ! 


-2.7-1. lx 
4.z+U.7a; 


«12 
6l2 


7.7-0.885x 
5.47+1.328X 


'r 25i ,r 25 „ 
™r 25 ,,r 25 „ 


-19.8-4.339x 
3.9-4.024X 


a r 2l ,r 2 „ 


0.3-1.511x 
-2-3.92(2! 


™ r 25' 


-7 379-2 148r 


C12 




n r 25i ,r 25 „ 


-0.112x 


9ri 2 ,r 2tl 


-10 5+5 5t 


I' 

1 15 


3.4+2.626x 


<2l2 


3.65+1.208x 


/ri",r 25U 


6+5.22x 


Sr 12 ,r 2; 


-4.5-0.854x 


m' r 

1 15 
<15 


-0.5+1.262i 
-10.392+0.2582: 


ar 2< 


-9+1.8191 
5-0.51x 


/r\i,r 25i 
/ri-,r 25i 


-5-2.666x 
-10-2.211 






l" 

1 25" 

m" 

1 25" 


-19-1.692i 
8+1. 119s 


ar i< 
ar iu 


10+4.171x 
0.5-0.992x 


/ri5,r 2 , 
/r 15 ,r 2 „ 


-19-3. 242x 
-2+21.9251 






n" 

1 25" 


-1.732+2.2131 
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TABLE V: Deformation Potentials. 

Si Exp. a EPM 6 LDA C k.p d 

b v -2.10±0.10 -2.12 -2.27 -2.27 

d v -4.85±0.15 -4.56 -4.36 -4.36 

(Sf + - a v ) 1.50±0.30 2.24 1.67 1.94 

+ i~£ - a v ) -1.4 -3.14 -3.03 

8.6±0.4 9 8.79 9.01 

S£ 15.9 13.85 15.1 

Ge Exp. a EPM 6 LDA C k.p d 

b v -2.86±0.15 -2.81 -2.9 -2.8 

d v -5.28±0.50 -5.31 -6 -5.5 

(H^ + - a„) 3.12 1.43 1.83 

(Sj + iSj - a v ) -2.0±0.5 -2.26 -2.86 -1.97 

9.91 10 10 

16.2±0.4 16.3 17 16.3 



"Cited by Ref. 6; b Chelikowsky and Cohen-based EPM 1 "' 53 ; c present work; d prcsent model. 
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TABLE VI: Eigenvalues and energy gaps at high-symmetry points in the Brillouin zone calculated 
with EPM, k.p and GW methods (see text for details). The state's SO splittings are shown in 
parenthesis. Averaged values over transitions between SO splitted bands are noted with a bar. All 



energies are in eV. 



ni „i „„ a 

states 


Exp. 


Si 
EPM 


G Wo 


k.p 


Exp. 


(je 
EPM 


G Wo 


k.p 


Li 


-6.8±0.2 6 ; -6.4*; 


-6.991 


-7.019 


-7.448 


-7.7±0.2 i ' 


-7.588 


-7.801 


-7.678 




-6.7±0.2 d 
















L 3 


-1.5 C ; -1.2±0.2 6 


-1.228 


-1.216 


-1.198 


-1.4±0.2 6 


-1.433 


-1.459 


-1.490 






(0.034) 


(0.033) 


(0.026) 




(0.187) 


(0.197) 


(0.188) 


Li 


2.06*; 2.1 s ; 


2.247 


2.095 


2.234 


0.744 d 


0.776 


0.64 


0.747 




2.4±0.15 9 
















L3 


3.9 ! ; 4.15±0.1 9 


4.324 


3.962 


4.245 


4.2±0.1 9 ; 4.4 C ; 


4.319 


4.227 


4.250 












4.3±0.2 d 












(0.016) 


(0.015) 


(0.007) 




(0.087) 


(0.103) 


(0.077) 


L 2 




7.334 


8.161 


8.031 


7.8 C ; 7.8±0.1 9 ; 


7.285 


7.495 


7.242 


Xi 




-7.711 


-7.823 


-8.087 


7.9 1 
-9.3±0.2 9 


-8.646 


-8.995 


-8.875 




-2.5±0.3 ft ; -2.9 9 ; 


-2.889 


-2.92 


-2.95 


-3.66 c ; -3.15±0.2 d ; 


-3.267 


-3.28 


-3.375 




-3.3±0.2 9 








-3.5±0.2 9 








^ 1 


1.10 , l.iU , 1.0 




1 221 


1.321 


1.3±0.2 d 


1.254 


1.045 


1.169 


Wi 


-8.1±0.3 6 


-7.512 


-7.653 


-7.662 


-8.7±0.3 6 


-8.512 


-8.88 


-8.638 






(0.006) 


(0.005) 


(0.295) 




(0.029) 


(0.025) 


(0.103) 


w 2 


-3.9±0.2 6 


-3.886 


-3.95 


-3.922 


-3.9±0.2 6 


-3.956 


-4.151 


-4.038 






(0.014) 


(0.013) 


(0.008) 




(0.042) 


(0.075) 


(0.154) 




-4.7±0.2 6 


-4.466 


-4.527 


-4.553 


-4.5±0.2 6 


-4.548 


-4.748 


-4.555 


Eg (A) 


1.17 d 


1.031 


1.076 


1.17 




1.04 


0.855 


0.961 


El(L) 


3.45 9 ; 3.46 e 


3.492 


3.311 


3.432 


2.05 d ; 2.22-f 


2.302 


2.099 


2.239 


E[(L) 


5.38 e ; 5.50 9 


5.577 


5.178 


5.443 


5.65 d 


5.889 


5.686 


5.750 


E 2 (X) 


4.32 e 


4.052 


4.141 


4.271 


AA5 f 


4.521 


4.325 


4.544 


"symbols 


of Rcf. 52; "as 


presented 


in Rcf. f 0; 


c as presented in Rcf. 42; 


d as presented in Rcf. 28; e Rcf. 43; 



^Ref. 44; 9 as presented in Ref. 45. 
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TABLE VII: Effective curvature masses and Luttinger parameters. 

Si Exp. G W EPM k. P d 

mf 0.9163° 0.925 0.89 0.928 

rrtf 0.1905° 0.189 0.198 0.192 

mf 1.8083 1.855 1.704 

m\ 0.1235 0.1535 0.131 

71 4.26°; 4.285° 4.54 4.01 4.338 
4.22°; 4.340 c 

72 0.38°; 0.339° 0.33 0.38 0.3468 

0.39°; 0.31 c 

73 1.56°; 1.446° 1.54 1.401 1.4451 

1.44°; 1.46 c 

Ge Exp. G W EPM k.p d 

mf 0.881 0.964 0.874 

rrtf 0.176 0.205 0.200 

mf 1.588^; 1.74° 1.626 1.763 1.59 

mf 0.08152' ; 0.079° 0.074 0.099 0.099 

71 13. O h ; 12.8 e ; 13.54 9.54 10.41 

13.25° 

72 4A h ; 4.0 e ; 4.32 2.75 3.045 

4.20° 

73 5.3 h ; 5.5 e ; 5.77 3.93 4.313 

5.56° 

"as presented in Rcf. 28; b l. Balslev and P. Lawaetz, as presented in Rcf. 22; c Rcf. 54; d prcscnt model; 
e Rcf. 55; f Ref. 56; 9 ReL 57; h Rcf. 17. 
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